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ABSTRACT 


The efficient use of a search-surface radar or sonar, in 
which one or more targetS appear on the screen intermittently 
usually demands a device for tracking the targets auto- 
matically. Such a device, called a "track while scan system", 
must make an estimate of each target's instantaneous position 
from the sampled-data information provided by the radar. 

For this purpose, an a-g filter and an optimal Kalman 
filter, that must track maneuvering targets, are analyzed 
here and compared in terms of tracking accuracy for tactical 


applications. 
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T” TNTEODUCTION 


Systems such as search-surface radar and sonar, in which 
the input data arrives intermittently, frequently require a 
device for continuously estimating the "present" value of 
the input. In radar terminology, this device is called a 
"track while scan system". 

In a more general sense, the term "track while scan" 
system, may denote any system which estimates the "present" 
value of a signal from the "past" sampled values of the signal, 
the sampling taking place at regular intervals. 

Vehicles, without maneuvering, of the class under consider- 
ation (such as aircraft, ships, and submarines) generally 
follow straight line constant velocity trajectories. If the 
vehicles were not able to deviate from these trajectories, 
1.e., could not maneuver, then the tracking problem could be 
solved quickly and simply using standard filtering such as 
ENEG-—S filter. 

Historically, a-f filters were designed to minimize the 
mean square error in filtered position and velocity under the 
assumption of small velocity changes between data samples. 
Thus, most a-ß filters have little capacity to track targets 
that either accelerate or maneuver (changing direction). An 
important early paper [1] defines 8 in terms of o with a a 
design parameter. This greatly simplified the use of maneuver 


detectors by reducing the design problem to a single variable. 
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A later advance in a-g filter design was achieved by 
minimizing the mean square errorin predicted position 
ai. 

A system that developed using the above technique is the 
VEGA LN I (Thomson-CSF, France). 

The final choice of parameters for the a-ß tracker will 
always end in a compromise between smoothing the input noise 
(measurement errors) and retaining some ability to follow a 
maneuvering target. Various performance measures have been 
used for this compromise. For example, steady state noise 
variance reduction and the ability to follow a target capable 
of responding to impulse accelerations are used as performance 
measures in [1] to derive good steady state filter parameters. 
Time varying noise variance reduction and the ability to follow 
a randomly maneuvering target are performance criteria which 
follow from [3]. In this paper, such a filter is analyzed, 
the frequency response, stability, noise characteristics and 
transient error are derived and plotted, and a scheme for 
Optimizing the two dynamic parameters (a-8) is suggested. 

A criterion for a-g filter design with a numerical and 
graphicai example is also given. 

Since the majority of tactical weapons systems requires 
that manned maneuverable vehicles, such as aircraft, ships, 
and submarines, be tracked accurately, an optimal Kalman filter 
has been derived for this purpose, based on the early work of 
[4,5]. The target model for tracking applications must be 


sufficiently simple to permit ready implementation in weapons 


ls 


mE 





systems for which computation time is at a premium yet 
sufficiently sophisticated to provide satisfactory tracking 
accuracy. 

The target acceleration selected for the Kalman filter, 
and hence the target maneuver is correlated in time; namely, 
if a target is accelerating at time t, it is likely to be 
accelerating at time t+t for sufficiently small 1. By sim- 
plifying the maneuvering model used in the Kalman filter above, 
the state vector can be reduced from six to four elements. 
The model simplification - simplified Kalman filter - is 
achieved by assuming (incorrectly) that the vehicle's change 
in velocity is uncorrelated between samples; i.e., the maneuver 
is white. In the comparison which follows, the Kalman filter 
and Simplified Kalman filter are compared with the o-8 filter 
in a variety of tactical environments with tracking sensors, 
utilizing Monte Carlo simulation techniques on realistic 


target trajectories to verify their theoretical performance. 
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II. FUNDAMENTALS OF TRACK-WHILE-SCAN 


Any system which performs a tracking function must obtain 
and utilize the basic target parameters of position and rates 
of motion. In the earlier system both position data and 
velocity data were used to maintain the tracking antenna on 
the target at all times, thus limiting the system to the classi- 
cal one target-at-a-time tracking function. 

In a track-while-scan system, target position must be 
extracted and velocities calculated for many targets without 
holding the radar antenna fixed on one target. Obviously in a 
system of this type, target data is not continuously available 
for each target track at a rate dependent upon the scan rate 
of the system. 

In a typical TWS system the data rate is one unit of data 
per second or a scan rate for the search radar antenna. 

Since the antenna is continuing to scan, some means of 
d ring and analyzing target data from one update to the next 
and beyond is necessary. The digital computer with its memory 
and computational capability is employed to perform this 
function and also: 

Gee TO provide a tactical picture 
A display of the tracks of all vehicles observed by sensors 
showing present positions, courses and speeds etc., is essen- 
tial for the deployment of a ship and its weapons. The computer 
enables rapid use to be made of sensor information thereby 


EMeuring that the picture is accurate and up-to-date. 
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Mn For In processes such as threat evaluation 
and weapon assignment 
Amer computer can forecast future likely positions of tracked 
vehicles and rapidly perform necessary calculations to assist 
operators in the assessment of threats and in the optimum 
weapon deployment to deal with them. 
(iii) To provide target information for weapon 

deployment 
The computer can be used to generate smooth tracks from noisy 
information thereby being able to pass information to a weapon 
sensor more accurately than information derived from a single 
plot on displav. Hence, the central concept underlying any 
TWS system is that the sensor itself continues to perform its 
primary function of search (Scanning) and data input while the 
remainder of the system performs the target tracking function. 
The sensor function simply provides target position data to 
the computer subsystem where target velocities and position 
prediction are calculated. In a military application the 
major advantage of a TWS system is the elimination of the 
process of target designation from a search radar to a fire 
SemerOl radar. 

The tracking information, developed in the TWS system, 
is used as a direct data input to the computation of a fire 
control solution. 

Therefore, as soon as a target is detected a fire control 
solution is available without the inherent delay caused by 


the designation process. The time required from first detection 
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to fire control solution is on the order of milli-seconds for 
a TWS ste as opposed to seconds or even minutes for a 
manually designated system employing separate search and fire 
control sensors. 

The focus of the following chapter has been to answer the 
question: "What functions should the TWS system perform in 
order to combine the search and tracking tasks into one 


integrated unit?” 


Pe A TRACK WHILE SCAN METHOD 
The method of solving the track while scan problem is 
based upon the assumption that the radar furnishes target 
position information once each scan. The scheme can be imple- 
mented by a combination of special radar circuits (hardware) 
and software. Existing systems also provide for operator 
modification of system tracking. 
Any track-while-scan system must provide for each of the 
following functions: 
(i) Target detection 
(ii) Generation of target acquisition and tracking 
"windows" or "gates" 
(iii) Target track initiation (assignemnt of targets 
to track files) 
(eee Target data input and track correlation 
LONE track window prediction, smoothing and 


posdtioning 
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l. Target Detection 
Target detection, localization and designation are 
accomplished by the radar sub-system in the usual sense. 
Hence, the TWS receives the following data from the outside, 
via wire links: 
(1) The surveillance antenna azimuth 
(11) The radar synchronization (presyne) 
(111) The surveillance radar video 
cem Generation of Target Acquisition and Tracking "Windows" 
A "window" can be defined as a small volume of space 
initially centered on a target, which will be monitored on 
each scan for the presence of target information. Each initial 
target detection will cause a relatively large acquisition 
"window" to be generated, centered on the position of the 
detection. When a target is initially detected the algorithm 
receives only the position data for that initial, instantaneous 
target position. The acquisition window is then generated, 


for example: 


Initial position = (Range, Bearing, Elevation) 
Range window = R + 1000 yars (914.1 meters) 
Bearing window = B+ 5° (0.0873 radians) 


Elevation window = E + 5° (0.0873 radians) 
WieMocquisieston window, Fig. la, is large in order to allow 
for target displacement during the following scans of the 
radar. If on the next radar scan the target is within the 


acquisition window, a tracking window is generated in the 
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Azimuth Window 
10 degrees 
0.1745 radians 


Range Window V 
2000 Yards A 
ee —d Elevation 
IT Window 
10 degrees 


0.1745 radians 








a. ACQUISITION WINDOW 


1.5 degrees 
— 0.0262 radians 






120 yards 
110 meters 








b. TRACKING WINDOW 


FIG. 1. TRACK WHILE SCAN VOLUMETRIC WINDOWS 
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same manner as the acquisition window. Although Fig. lb 
shows only the very small, 120 yards, 1.5°, 1.5? tracking 
window, in actual practice intermediate window sizes are 
generated until a smooth track is achieved. 

A rack: Initiation 

Concurrent with the generation of the acquisition 
window a track file is generated in order to store the posi- 
tion and window data for each track. In addition to the basic 
position and window data, calculated target velocities and 
Beeelerations are also stored in each track file. Track files 
Ame stored within the digital computer (or processor) sub- 
systems memory and the data is used to perform the various 
calculations necessary to maintain the track. 

Each track file occupies a discrete position of the 
digital computer's (or processor's) high speed memory. As 
data are needed for computation or new data are to be stored, 
the portion of memory which is allocated for the required data 
will be accessed by the system programs (software). In this 
manner a diversity of data in addition to the tracking data 
may be stored in the "track" file, for example, ESM data, IFF 
information. The generation of the track file begins with the 
initial storage of position data along with a code to indicate 
that an acquisition window has been established. If target 
position data is obtained on subsequent scans of the radar, the 
file is updated with the coordinates; the velocities and 
accelerations are computed and stored, and the acquisition 


Window code is canceled. 
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The acquisition window is then decreased in size to 
the tracking window and the track code is stored which indi- 
m S an active track file. 

As the radar continues to scan, each input of data 
is compared with the window positions of active track files 
until the proper file is found and updated. However, it should 
be noted that the search for the proper track file is generally 
not a sequential one-to-one comparison. This method is much 
too slow to be used in a system where speed of operation is 
one of the primary goals. 

This idea of comparing output data with window posi- 
tions leads us to the problem of correlation of data input to 
e files and what to do if correlation is ambiguous. 

fee resolution of Track Ambiguity 

Track ambiguity arises when either multiple targets 
appear within a single track window or two or more windows 
Overlap on a single target. 

This occurrence can cause the system to generate 
erroneous tracking data and ultimately lose the ability to 
maintain a meaningful track. If the system is designed so 
that an operator initiates the track and monitors its progress, 
the solution is simply for the operator to cancel the erroneous 
track and initiate a new one. 

For systems which are to be automatic, software (pro- 
gramming) decision rules must be established which will enable 


Ee tracking program to maintain accurate track files. Decision 
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rules as to target input data and track correlation can be 
logically developed with an understanding of the definition 
of the track ambiguity problem. 
5. Track Window Prediction, smoothing, and Positioning 

In a track-while-scan system tracking errors also 
exist due to target motion. The tracking window now has 
replaced the "tracking antenna" and this window must be posi- 
tioned dynamically on the target in a similar manner as was 
the "tracking antenna". However, there is no "servo" system 
tO reposition and smooth the tracking window's motion. This 
repositioning and smoothing must be done mathematically within 
the TWS algorithm. To this end, smoothing and prediction equa- 
tions (Eqs. 2.1-2.3) are employed to calculate the changing 
position of the tracking window. Instead of the system "lagging" 
the target the tracking window is made to "lead" the target 
and smoothing is accomplished by comparing predicted parameters 
with observed (measured) parameters and making adjustments 
based upon the errors derived from this comparison. 

Figs. 2 a-b and 3, illustrate the TWS principle, 


track positioning and the TWS general processing loop of the 


alpha-beta type, respectively. 


B. DEFINITION OF TRACKING EQUATIONS FOR a-B FILTER 
Tracking in a track-while-scan (TWS) system consists 
basically of two functions: 
a. Smoothing 


D correlation 
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DIGITAL COMPUTER OR 
PROCESSOR 
(Software) 





a. IWS PRINCIPLE 


ACTUAL TRACKING 
APPROXIMATION 


POSITION 
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b. TWS POSITION TRACKING 


FIG. 2. TWS PRINCIPLE AND POSITION TRACKING 
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Smoothing is the processing of the sensor reports to derive 
an estimate of both target velocity and position. 

Correlation is the sorting of sensor reports into groups, 
to determine which belong to which target then predicting a 
new coordinate on which to center the correlation region. 

The constants of proportionality, a and g, used in 
correcting the position and velocity, respectively, of the 
estimated target course, completely characterize the perform- 
ance of the TWS system. These constants are the dynamic 
parameters or so-called "Smoothing constants" of the system. 

Finally, the simplest case target tracks are based on 
smoothing and prediction of an alpha (a) - beta (8) tracker 
Operating in a cartesian coordinate reference frame. The 
information is ordinarily obtained from a coordinate converter 
SeeLacing on the raw polar to cartesian transform. The o-8 


filter described by the well-known a-8 tracking equations [6,7]: 


N N N N : : 
Xo = XO +t a (Xy - Xp) smoothundgeEaguatrons 2 
N N 
ne IS 
Ve = Vo + nd prediction Equation (2.2) 
A = SE: + a T prediction Equation 2") 
where: 
Xo = smoothed position 
Xp = predicted position 
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V = smoothed velocity 


S 
Zu = measured position 
a,B = filter parameters or filter gains 
T = sampling time or time between detections 


EnbSEPeuting Eq. (2.3) into Eqs. (2.1) and (2.2) yields 


No Ca qN-1 I N-1 N 
Xc = (il a) Xo + (l-a) TVo + ax, (2.4) 
N | | 8 ,N-1 ERN 2 8 sN 
Vo = T Xo + (1 8) Ve T T Xy (2:5) 
N+1 N 
= 2.6 
X. Xg + Vs T ( ) 
Or 
" 8 o N-1 
Ko (l-a) (1-a)T Xç 
= + 
V x (1-8) V 
ms! T S 
a fès lè 
+ (2 7) 
Ë 
T 
and 
N+1 ZEN 
[Xp] = [1 T] Xo (2.8) 
Vo | 
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These equations represented in block diagram form appear in 


Big. 4. 





Fig. 4. Block Diagram of an a-8 Tracker 


where the following matrices and vectors are defined: 


WE MES a 
aS Bx 
-& (1-8) B/T 
E A 
C = [1 T) x = 
ioni 


1. Stability Analysis 
Since the system characteristic equations are in the 
form of difference equations, Z-transform theory is helpful 


in determining filter transfer functions and stability criteria. 
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For the system of the form: 


Xo] - AR + B (Xm = CX) 


Or 


>< 
Il 


N+1 a BC) Xy MB 


thus, 


= oq (A-BC)] | B (2.9) 





Meang Eq. (2.9), the Z-transform for each of the several 
components of Eq. (2.7) and (2.8) has been computed and is 


given below. Hence, the a-ß tracker transfer functions: 


x (2) 
S ala (zel) Eg) 
GAZ) = = ——N.o“ak . + -— Smoothed (2.10) 
“e Xy (2) D(2) position 
B 
V (2) pi (2-1) 
Gaz) ^ — = ——n Smoothed (2.11) 
vs Au EL SEH velocity 
X (Z) 
P a (2-1) + 82 . 
Gn (Z) = pia Predicted 2.12) 
ns XM 2 CU, position 
B 
VE (2) p (2-1) 
G.(Z) = == = —— H Predicted (2 12) 
ME Xy (2) EIS, velocity 
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where the characteristic equation of the system is 


Dia) ZA) =. 6 


or 


DI (2-54 (l-a) = O (2514) 


The stability of the system can be determined from the loca- 
tion of the roots of the characteristic equation above. 

A few methods are available for determining whether 
or not a polynomial in Z contains a root or roots on or out- 
side the unit circle. One method is to modify the Routh- 
Hurwitz stability criterion. The Routh stability criterion 
tells whether or not any of the roots of a polynomial lie in 
the right half of the complex plane. 


Since the following transformation 





Dapssche interior of the unit circle in the Z plane to the 
lef-half r plane, with this transformation, permits the 
application of the Routh stability criterion to the polynomial 
in r as in continous-time systems. 

Now Eq. (2.14) becomes a polynomial in r, in the 


form 
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gri pne 4-2; -B = 0 (2215) 
So, the resulting requirements for stability are 
Cem; pavo and (2a + B) < 4 2,156) 


An additional stable condition exists when g = 0. Thus the 
resulting necessary and sufficient conditions for the stability 


of the track-while-scan system are 
BCEE ME and (2a + 8) < 4 (2717) 


These inequalities determine a "Stability triangle" in the 
a-8 plane, for which all internal points and all points on 
the base (g = 0) in the interval 0 < a < 2 correspond to a 
Stable system. This triangle is shown in Fig. 5. 

The conditions for underdamped, critically damped, and 
overdamped transient response are found by inspecting the 
men cr Ene discriminant of Eq. (2.14). The resulting condi- 


tions are: 


EE CSS 
PS 1, EE = 48 > critically damped 12519) 
a « 1, B < 1, (a+8)° > 48 <= overdamped (2.20) 
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B) 


( 


BETA 





ASYMPTOTICALLY STABLE 
2a +8 < 4 


21 





All other values of The transient response 


(0,8) inside the stability contains at least one 
Errangle-and on thé base <=  demped oscillatory natural 
(8 = 0) in the interval mode with a rate of 

DI 2 oscillation equal to one 


half the sampling frequency 


The permissible values of a and 8 are shown in Fig. 6 


2 


1 2 
a 


Fig. 6. Allowable Values of a and $8 in a-8 Tracker 


Further restrictions will be placed on this region, when the 
frequency response characteristics of the system are examined. 
D new Response Characteristics 
In this section, the frequency response characteristics 
of the a-g filter are formed. This approach has not been 
done before, in detail. 


The frequency response of the filter can be found by 


placing Z = Li and pen COS Sin we into Eqs. (2.10)- 
(2.13). These equations are now in the form: 
a. JE E [a (cos2uT-cosuT)+8coswT] + ] [a (sin2uT-sinuT)+ssinuT] (2.21) 


mee) 


52 





sur (5 (cos2uT-cosuT) jose (É(sin2uT-sinuT) ] 


(e- `) - (2522) 
Sus Die): 
E E > [(a+8) cosuT - a. j [(a+8) sinuT] (2.23) 
D (el ) 
| [=(cosuT-1)] + 3 sinn 
jwi T Joe 
Gp (e ) eee VO — (2.24) 
D(e- `) 
where 
Deu, =  [cos2wT- (2-a-8)cosuT+(1-a)] + jisin2uT=(2=a=8) sinurT] 
The amplitude and phase characteristics of G (smoothed 


XS 
position) and Gyp (predicted position) are plotted and shown 


in Figs. 7-16, for several values of a and 8. Also the ampli- 


tude and phase characteristics of G (smoothed velocity) and 


VS 
Gyp (predicted velocity) are shown in Figs. 17-26. All the 
amplitudes and phases are plotted as a function of a, 8 and 
WES 

Observing these figures and Eqs. (2.7) ~ (2.8) one 
fines that us is the result of passing Yen through a low pass 


filter, Ve is the result of differentiating e and a should 


M’ 
never be larger than one. 

It may be seen also, for positional smoothing, if 
a = 0, all sensor information is ignored whereas if a = 1, 
there is no smoothing of positional information. Similarly 


8 = 0 causes sensor information to be ignored in the estima- 


tion of velocity whereas g > 1 will cause overcorrection. 
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Considering the figures for several values of a < 1, 
a seems to control the bandwidth of the low-pass filter and 
8 has more control over the damping. In fact g should be 
somewhat smaller than a such that resonant spikes do not 
Becur. 

Better accuracy occurs between smoothed and predicted 
positions, when a has very large values compared with 8. 

Table I, which follows, summarizes the accuracy and 
provides an important design tool. In general the frequency 
of the input signal, the sampling time T, and the filter 


parameters o and 8 control the filter's response. 


54 











ACCURACY (0-1) BETWEEN SMOOTHED AND PREDICTED 
POSITION FOR wT UR o um = ur 
CL 

x 0.1 0.5 0.75 0.9 
Oro 1.0-0.90 1.0-0.98 1.0-0.99 
0.05 1.0-0.60 1.0-0.90 1.0-0.94 
DO 1.0-0.33 1.0-0.82 1.0-0.89 
0.25 100.71 
0.50 1.0-0.50 
0.75 | 1.0-0.33 


ACCURACY (0-1) 





BETWEEN SMOOTHED AND PREDICTED 


VELOCITY FOR wT = 0 to uT = T 
A x DET 0.5 0.75 0.9 

poll 1.0* 1.0* IO 

0.05 o 1.0 120 

0.10 1.0 rO 1.0 

0.25 7 

0.50 1.0 

0.75! 4: (6 


“As we expect 100% accuracy, 


since the vehicle was 


generally assumed to follow straight line constant 
velocity, when an a-ß filter is used. 


ACCURACY BETWEEN SMOOTHED AND PREDICTED, 
POSITION AND VELOCITY 


TABLE I. 
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Wie OSE CHARACTERISTICS OF a-8 FILTER-ERRORS-CRITERIA 


In the study of any filter it is essential to know the 
characteristics of the desired signals and the noise which 
excites the filter. It is desirable also to know how the 
choice of a and 8 affects the degree to which the noise 1s 
smoothed or exaggerated by the system. The description of the 
noise processes, prediction errors and methods (criteria) 


for designing an a-8 filter, proceeds as follows: 


FR NORMANIZED NOISE POWER OF PREDICTED POSITION 


By making the following change of variable: 


D OR E TE ES an 


ESNeIons4(2.7) and (2.8) describing the filter can be 


rewritten as: 


N R e N-1 _ N 
Xo (1-0) (l-a) Xe a [Xy] 
E | + (3 2) 
aX | -8 SS | 8 | 
Nel | TE 
[Xp] = ji 1] Xo i 5 
I-A X 


where now 
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The mean and 


where 


14 


12 


22 


AM 





M idis a | 
| a. 
| 28 ER |B | 
el 
EO Ze 
AX 


atx i B'Xy (3.4) 
a'e ta'l 4 rta (3.5) 

Rèn , the expected value, 

E o 

Es 2. 


E[(AX - AX) (AX - AX)] 


Standard deviation of measurement error. 
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The covariance equation for the filter becomes 


- - N - 
Pry x E 2 (1-0) 
Pi = -8(l-a) IB San MIA AO 
| 
| 
2 2 

E | 8 A) GE 

a“ C 

a aB | 
SCH 


(3-6) 


The variance of the predicted position is easily computed 


in terms of the variances from Eq. (3.6) 


pie) * ze, Hin 


The steady state solution of Eq. (3.6) is computed by leting 


P(N) equal P(N-1) and the resulting algebraic equations are 


solved; 





A AS 
2 u a(4- 2a = B) 
XM 

Mo Mo) 

E “ oa(4 - 2a - 8) 
XM 


DIS 


(38) 


29) 








P 2 3 








22 Bän eo + 28 - a8) 
ER E a(4 - 2a - 8) Geck 
XM 
The predicted noise power is normalized 
52 
R = 0° = — (Sail) 
7XM 
iene inglesi. (3.7) (3.10) into Eq. (3.11) yields 
52 2 
s XP _ a (2 28 -a8) + B(29 + a - ag) 
E oos. "7 so. c pesa: 
IXM 


Hence, the predicted noise power is plotted as a function of 
a and 8 in Fig. 27. This figure shows that by appropriately 
adjusting o and 8 the input noise power can be reduced. One 
also finds that the parameter 8 greatly affects the predicted 
noise power, as expected, since differentiated noise is quite 
"noisy" and 8 affects this quantity. 

An interesting phenomenon occurs for small values of 
a and large values of 8. The noise power increases sharply. 
The reason for this may most easily be seen by looking at the 
frequency responses in Figs. 7-16. The large resonant peaks 


in the response allow a lot of noise to come through. 


B. PREDICTION ERROR DUE TO MEASUREMENT ADDITIVE NOISE ERRORS 


In the a-8 tracker the observations are in the form 


_ N 
x = Xu + v ESB) 


DE 
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FIG. 27. PREDICTED POSITION NOISE POWER AS A FUNCTION 
OF a AND 8 
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w 





where: 





= additive measurement noise 
We - M 
= V y IR 
= 0 LO N # M 
and by letting 
Variance p = VAR (XP) 
Variance Reduction Factor = VRF 
hence 
N+1 VAR (Xp) 
VRF (X ) = ro (3m 
P 2 
G 
` 
OI 
N+1 
VAR(X ) 2 
ee eB e TE) 
D 2 a (4-2a-8) 
g State) 


mne VRE is plotted in Fig. 28 as a function of o and 8. Com- 
paring this figure with Fig. 27 (normalized noise power of 
predicted position), the two figures look alike, but this is 


not true, when both graphs are plotted together in Fig. 29. 


See eee CTO ERROR DUE TO CONSTANT ACCELERATION 


Constant acceleration X(t) results in fixed error in 


PL 


predicted position Xp 


on reaching steady state, i.e., for 


Die oos 
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NORMALIZED PREDICTION ERROR DUE TO ADDITIVE NOISE, VRF 
(VARIANCE REDUCTION FACTOR) 


PIG. 28. 
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PREDICTION ERROR DUE TO ADDITIVE NOISE AS A 
FUNCTION OF a AND 8 
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FIG. 29. PREDICTED POSITION NOISE (R) AND PREDICTION 
ERROR DUE TO ADDITIVE NOISE (VRF) AS A 
FUNCTION OF a AND B 
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Specifically for noiseless constant acceleration trajec- 


tories, 
X X (3.16) 


the steady state prediction error becomes 


. 2 
s EI Gs, 


N 
€ = lim (Xy SE TI 


N^ o 


* 
Common names for steady state prediction error e due to 


constant acceleration are: 
(1) Dynamic error 
Cle Truhcatron error 
(111) Systematic error 
(iv) Bias error 
Typical design procedure for balanced dynamic and random 


errors is 


3 (Standard deviation, E (Dynamic, 
of random error error 


Chat 15 


3 | var (x+) Ben, (3.18) 


D. TRANSIENT ERROR 
The transient error of the system becomes significant when 


the input is switched from one time-shared target to another 
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or when the target makes a sharp maneuver. Since o-8 trackers 
are initialized by original estimates of position and velocity 
and in many applications these estimates are so inaccurate that 
the filter exhibits transient errors before settling to a 
steady-state solution. In the transient phase, when the esti- 
mates are poor, one would like to rely on the measurements 

and thus weight the errors more heavily. After several measure- 
ments have been made, the accuracy in prediction increases. 


So, the transient error yields 


2 
Y N+1 N 2 Be T' (2-0) NOS 


N+1 D M ~ aB(4 - 2a - 8) 


The transient error D is plotted and shown in Figs. 30-32, 


N+1 
Boma function of o, B El T. Observing these figures, one 
finds, that the sampling time T greatly affects the transient 
error. For large values of T the transient error increases. 


An interesting phenomenon occurs for small values of a and 


small values of 8. The transient error increases sharply. 


BERO :MALTSTEADY STATE RELATIONS FOR THE a-8 TRACKER 

The simplest possible system is one in which a and 8 
are fixed constants. The theory of such systems has been 
studied in [8]. For such systems, T.R. Benedict and G.W. 
Bordner [l] applied the calculus of variations to derive the 


following relationship: 
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B = (3720) 





A few years later S.R. Neal [7] used the results of linear 


estimation theory to derive the relation: 


(de LS = 88 (3:21) 


As a matter of interest, the above Egs. are plotted and 
compared in Fig. 33. Observing these figures it is inter- 
esting to note that the relations (3.20) and (3.21) are quite 
n i or for 0 « «e «0.4 and 0 < pg < 0.1. 

Constant parameter systems suffer from the incompatible 


demands that good smoothing requires heavy damping (i.e., small 


values of a and 8 - small noise response, Figs. 27-28), while 
good response to maneuvers requires light damping (i.e., large 
values of a and g - small transient error, Figs. 30-32). Light 


damping and therefore large noise response, can lead to low 
probabilities of weapon sensor acquisition and problems of 
plot-to-track association. Heavy damping, implying poor 
maneuver response (large transient error), can cause sudden 
loss of tracks (sometimes termed track death) through failure 
to associate with subsequent plots. 

These limitations led some workers to optimize for varia- 
ble parameter systems, where a and 8 are varied according to 
the state of the track. Some systems have been developed 


wherein a and 8 were initially selected arbitrarily and changed 
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during program development by trial and error. Various 
Operational sets of values being derived for various states 
of track. Such methods are adaptive and are usually economi- 
cal in computer use, both in terms of required storage and 


run time, but generally have no theoretically optimum adaptation. 


fae DESCOUNTED LEAST-SQUARES CRITERION 


To find linear trajectory 
Xp = X +NT Vo (3722) 


which minimizes sum of weighted errors between xh and 


2 
Ay Ay , Zu ee. 


e. = ] DN NE 3223) 


Ec 201 str too many degrees of freedom for 


selecting gain terms, a and g. Minimize 
kai N_,N,2,0 N-=1_,N-1,2,1 N-2 .N-2,,2 
(ey? = (Xy x) 0° + (Xy Xp ) e” + (X Xp Ja” + 
(3.24) 


by solving the above equation in order to weight the differ- 
ence between measured and predicted values. This yields 


in simple gain terms 


(== o) (325) 


e 
Il 

















Ee — S io = O) (3.26) 


Common names for discounted least-squares a-8 filter are 

(i) Critically damped a-ß filter 

(ii) Fading-memory polynomial filter of degree 1 
Equations (3.25)-(3.26) are plotted as a function of (theta) 
mm Fig. 34. 

More recently, processes have been developed in which 

and are made to change with time in order to continually 
compute the least-squares line through the observations. Such 
approaches assumed that errors are equally distributed in 
x and y and had a constant standard deviation. The formulae 
Re changing a and g in this manner were worked out in [9]. 


Por the incorporation of the nth measurement: 


Ac 1) 


Q IGES n (3.27) 
dl n (3.28) 
N(N + 1) ° 


The above equations are plotted as a function of the number 
of measurements N, in Fig. 35. 

It is clear that, for large N, a and g tend to zero, 
i.e., observations will be increasingly ignored. This sug- 
gests that there should be some maximum value of N. The 
maximum value used is generally 7 to 15. Such a method may 


be made adaptive if a means of detecting changes in motion 
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is used, i.e., if turn detection is provided. Then, if a 
turn is detected, the values of a and 8 may be raised simply 
by lowering N. Doing this will improve the turn following 
capability. 

1. Numerical Example 


Assume the following; 


Range measurement = CR = SIOE 
Eer = 5g = 150 ft/sec 
Desire | VAR GS) ES 


etici table eriticallv damped a-8 filter design. 


Brissubsertucıing EGS. (3.25) and (3.26) 


a = l- gè 
2 
B = (l - 38) 
ERE. (3.15) 
N+1 
n P M 2a. * as * 26 
zè a(4 - 2a - B) 
R 
and letting 
VAR (xy) rc oz ES 


75 





yields 


VRF = 1.074 

8 = 0.5 
hence 

Ques 0.70 

poe 0525 


Brom Eq. (3.18) 


* 
a VAR (XN) EO = oa e ft 


Ste from Eq. (3.17) 





hence 


T = 0.3848 sec. 


Fema liv, From Eg. (3.19) 


= Tê (2 - a) 
N+1 ag(4 - 2a - B) 
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yields 
= E 2 


The values for a = 0.75 and 8 = 0.25 satisfy the critically 


damped Eq. (2.19) 
2 
CHEN) = 48 and B < 1 


The above values for a-B, VRF and D can also be determined 


bu using Figs. 27-31-34. 
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iy eo ead LNG OPTIMAL TRACKING FILTER 
FOR MANEUVERING TARGET 

The work which follows is an attempt to apply optimum 
filter theory to the tactical radar environment. There it is 
desired to obtain optimal estimates of target positions and 
their predicted tracks. 

There are several possibilities for structuring a digital 
filter. From the work of Kalman [3], the filter should con- 
tain the same dynamic model as that of the incoming signal 


shown in Fig. 36. 
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Fig. 36. Schematic of Filter Structure in Relation to Signai Process 
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The modeling of the various system components involved 
in a tactical weapon system, such as the radar meaSurement 
process and the target itself, is essential in the design 
of practical tracking and control algorithm. By modeling 
the target to be tracked and the accuracy of the radar's 
measurements [10], then a practical tracking procedure, con- 
Sistent with the computer limitations and weapon system 


requirements, can be designed. 


A. SENSOR AND VEHICLE MODELING - DYNAMIC EQUATIONS 

The tracking systems under consideration utilize sensors 
that provide measurements of range and bearing. The selection 
is intended to reflect that this pair of measurements is most 
common, however, other output measurements as elevation and 
range rate (Doppler) are often available. 


The vehicles to be tracked can be modeled by the state 


equations: 
X(K+1) = oX(K) + Gu(K) (4.1) 
where 
Y (K) range at time K 
XE) range rate at time K 
E 6 (K) Ë bearing at time K 
8 (K) bearing rate at time K 


-— 


= vehicle state vector at time t - KT 
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u. (K) change in vehicle range rate between 
time K and time K+1 


u(K) = = 
u, (K) change in vehicle bearing rate between 
time K and time K+1 
1 T 0 0 
0 d 0 0 
e — ] 1 1 A 
$ 0 0 1 T state transition matrix (o) 
0 0 0 T 
Mor W 
G = 1 0 
0 0 
0 Hi 
and 
T = sampling period. 


The tracking sensor measures target position, range and bearing 


or elevation and provides the following output equation: 


measured rate at time K 


Y(K) = HX(K) + V(K) = 
measured bearing at time K 
(4.3) 
where 
1 0 0 0 V4 (K) 
H = and V(K) = (4.4) 
0 0 1 0 V, (K) 
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The measurement noise covariance matrix R(K), satisfies 


a 0 
A R (K) 
R(K) = E[V(K)V (K)] = (4.5) 
0 me 
8 (K) 


assuming the noise V4 (K) and V, (K) are independent. The 
selection of sensor coordinates (R,8), at this point, for 
the covariance matrix R(K) has been made because the output 
matrix H assumes the extremely simple form shown, and the R(K) 
becomes diagonal. When Doppler measurements are avallable, 
this selection of sensor coordinates becomes extremely advan- 
tageous because the cartesian forms for H and R(K) becomes 
complex and time varying and often impose computational 
penalties for real-time implementation. 

Consider the following polar to cartesian transformation 


me., from R,9 to X,Y) 


XxX = R sin 9 


(4.5a) 


K 
Il 


R cos ° 


For the Kalman filter in XY coordinates, the measurement 


covariance matrix R(K) is a function Of the radar-target 


geometry. Therefore, the elements of the covariance matrix 
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XX 
R(K) = 
2 
“xy 
are 
2 u 2 2 2 
xx = CR cos 808 + R 
2 2.2 2 
Oyy = SR sin 8 + R 
2 = 2 2 
Oyy = lo, - R 
2 2 
where 7n and d 
bearing 


are the variances of the range 


(8) measurement errors, respectively. 


„2 
XY 
(4.5b) 
po 
YY 
2 P2 
Ca sin 8 
0% SETE (405) 


2 
Cal Sin 0 eos 9 


(R) and 


It is inter- 


esting to note that in general, the coordinates after the 


transformation (X and Y) are not independent. 


cases where they are independent occur for 6 = 


G 
or for R = —. 
Mo 


O e. 


The singular 


090, 


These special cases are easily 


explained from the assumption that the measurement errors are 


Gaussian, with typical contours of equal probability, which 


are ellipses. 
Hence, 


the frame of reference, 


are zero and the independence of X,Y errors for 8 = 0°, 


etc. 


When R = 


the covariance terms R(K) 


if the axes lie in the directions of the axes of 


R(K) 


E 2l 


909; 


OR/ Ca these ellipses reduce to circles which 


may be considered as limiting ellipses whose major and minor 
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axes lie along the axes of the reference frame. Finally, 
since the observation y(K) consists of x and y position 
measurements, hence, the observation V(K) can be assumed to 


have covariance 


gè 0 
= R 
R(K) = 2 (4.5d) 
0 Op 


B, STATISTICAL DESCRIPTION OF TARGET MANEUVER 

The input sequence u(K) is additive state (or maneuver) 
noise that results in the vehicle deviating from a constant 
velocity trajectory. 

Although the maneuver history and observation noise are 
not independent, the covariance Elu(K)v” (3) 1 is zero. Indeed, 
the radar cross section of a piloted vehicle changes during 
a maneuver, causing the radar observation noise to depend on 
the particular maneuver being exercised by the vehicle. 

The maneuver noise is neither Gaussian nor white. For 
example, the pilot of an aircraft moving at constant velocity 
will generally not maneuver unless threatened by either radar 
detection or attacking vehicles. His maneuver will then often 
be a turn or an increase or decrease in his forward velocity. 

A typical maneuver probability density is shown in 


ERRO 37. 
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Fig. 37. Typical Probability Density of Target Maneuver 


The quantity A is the maximum acceleration which the plane- 
pilot combination can withstand. Values of the density 
between non-maneuver (u = 0) and maximum (u = +A) are non 
zero because: 
(1) The vehicle may not be accelerating at its 
maximum rate 
(11) The projection of a circular maneuver on any 
dimension can give values of u from -A to A. 


Clearly, then, the maneuver density is not Gaussian. 


CRD TSERETE TIME EQUATIONS OF MOTION 

It is often desirable to whiten the maneuver noise so 
that system equations to which optimal filtering theory 
applies can be obtained. This is done, as in [4]. The 
whitening procedure for a discrete signal is analogous to 
the procedure developed by Wiener and Kolmogorov to white 
continuous signals [11]. So, u(k) may be expressed recursively 


in terms of the white noise sequence W(K) by 
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u(K+1) = pu(K) + W(K) 


(4.6) 


Using the above equation, the following set of system equations 


are obtained having white noise sequences W. (K) and W, (K) 


as their only inputs: 
X(K+1) = 4X(K) + GW(K) 
Y(K) = HX(K) + V(K) 


where 


r (K) ] 
r (K) 

u, (K). 
6 (K) | 
A (K) | 
u, (K) 


and 


W(K) = 


85 


(4.7) 


(4.8) 


iio m) 


(4.10) 





Also 


1 T 0 0 0 0 
0 1 1 0 0 0 _ 
0 0 o 0 0 0 State Transition 
Dues m 
H 0 l : 0 Matrix 
0 0 0 0 1 L 
ON NN 0M (4.11) 
0 0 
0 0 
1 0 
e = 0 0 Cae 2) 
o 0 
lo 1] 
1 0 0 0 0 0 
H = 
0 0 0 ] 0 0 | 
= [Observation matrix] (4.13) 
and where 
NN. 0 
T M1 
Q(K) = E[W(K)W (K)] = 
2 2 
0 Oy (1-9 ) 
Covariance of 
= the measured 117) 


noise 
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! 2 . 
The values for the maneuver variances mi and the correlation 


coefficient, p, depend on the maneuver characteristics of 
the vehicles being tracked. 

Hence, the variances, o Of target acceleration are 
calculated using the model illustrated in Fig. 37. The 
target can accelerate at a maximum rate A (-A) and will do 
1 and a probability P. Or not 


accelerating at all, with an assumed uniform probability 


each with a probability P 


distribution of amplitude: 


l - (2P, +P.) 
Eee cm 2 (4.15) 


2A 
ef accelerating between -A and +A. 
The acceleration variable, therefore, has mean zero and 


variance 


sf, = É (4.16) 


Consequently, the variables u; and u_, which are assumed 


D 


independent, have zero means and Ui has variance 


ema + 4P, - P) 
c? = EL 1 2. (427) 
Ml 3 
while u, has variance 
AD“ (1 + 4P, - P.) 
zè e ME SEN (4.18) 
M2 2 3 
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where R is the target range from the sensor and where all 
these quantities have appropriate units. 


The correlation coefficient (|p| < 1) can be modelled 


by 


l- uT r " 
= EREDI _ T 
SE 
Mi 0 T> Ł 
u 


The quantity, u, is essentially the inverse of the average 


Maneuver duration. This correlation model is analogous to 


92 e aT 
Mi 4 


where a is the inverse of the continuous maneuver time con- 


that in [4], which has the discrete time form r(K) = 


stant of the target. Hence, p equals Pa When aT is small, 
p can be approximated closely by l-aT, so that a and uy become 
identical. 

The two extreme cases occur when o is unity and when p 
approaches zero. So, the first case represents the completely 
correlated case, and the second is the completely uncorrelated 


case, namely white noise. 


ae Lee DESCRIPTION 

Three types of filters are considered as potential can- 
didate algorithms for tracking vehicles that are described 
by the model just discussed. These filters are the Kalman, 


the simplified Kalman and an a-ß filter. 
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l. The Kalman Filter 
The Kalman filter uses the augmented version of the 
model presented earlier in order to obtain white excitation 
(maneuver) noise. 
The method of computing the optimum estimate (the filter) 


is as follows: 


SEET (4.20) 


X(K|K) = X(K|K-1) + P(K|IK-1)HT(HP(K|K-1)HT+R(K)] À 


(fan 2 
[y (K) -HX (K| K-1)] 


where 
T T 
P(K|K-1) = éèP(K-1|K-1)d +G0(K-1)G (4.22) 
T T -1 
P(K|K) = P(K|K-1)-P(K|K-1)H (HP (K|K-1)H +R(K)] “HP(KIK-1) 


(4:23) 


In these equations 


(i) The double argument always denotes an estimate 


(ii) X(K|K) = Filtered estiamte of X(K) 
(111) X(K|K-1) = one-step predicted estimate of X(K) 
(iv) P(K|K) = Covariance of filtered estimate 
(v) P(K|K-1) = Covariance of predicted estimate 
(vi) Ly (K) -HX (K|K-1) ] = Error in the predicted 
observations 
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H 


Il 


(vii) G(K) P(K|K-1)HT [HP (K|K-L)HT+R(K)] 


= Kalman Gain Matrix 


= Is a matrix of adjustment coefficients. 
The matrix, G(K), reflects the relative 
confidence one should have in the 
observed data as compared to the 
predicted estimate. 


The filter iS initialized on the basis of two observations 


as follows: 


y, (2) 


d 


x (212) E (1551) 


yo (2) 


1 


The elements of the corresponding covariance matrix, P(2|2), 


aros 
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2 2 
7n geil 0 0 0 0 
2 
20 
2 2 R 2 
oy T Cmt e POMI 0 0 0 
2 2 
0 POMI 9M1 0 0 0 
P(2|2) = 
0 0 0 0% Be, 0 
0 0 0 ET we 2% 
M “M2 z2 iis pP 
2 2 
| 0 0 0 0 PSM) Ou: 
(4.25) 


2 | 
M17’ °M2 should be calculated as discussed 


previously (Eq. mi An a 


The quantities ae 


A block diagram of the discrete Kalman filter is shown in 


big. 38. 


Yy 





Beie © Xelg-1 * eU O 
sss ESI 


Fig. 38. Block Diagram and Equations of the Discrete Kalman Filter 
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2. The Simplified Kalman Filter 

By simplifying the maneuver model used in the Kalman 
filter, the state vector can be reduced from six to four 
elements, and the number of independent components of the 
covariance matrix from ten to six. The model simplification 
1s achieved by assuming (incorrectly) that the vehicle's 
change in velocity is uncorrelated between samples, i.e., 
the maneuver is white. 

The regular Kalman filter requires two augmented state 
variables in order to whiten the target maneuver. If the 
maneuver iS assumed white, no augmentation need be performed 
and the simplification just discussed occurs. This simpli- 
fied Kalman filter can therefore also be referred to, as an 
unaugmented Kalman filter. 

The utility of this filter is greatest, therefore, 
when either sensitivity of tracking performance to assumed 
Maneuver correlation is small, or when the target maneuver 
approaches whiteness relative to the sensor data rate. 

The equations for this filter and all quantities 


except the following have previously been defined: 


2 
“Ml U | S 
Covariance of 
Q(K) = E "assumed white" (4.26) 
2 maneuver noise 
0 M 


and the elements of the corresponding covariance matrix, 


P(2|2), are: 


DIO 











2 
On op/T 0 0 
. 2 
20 
2 2 R 
0/1 OM] * 22 0 0 
P(2|2) E 
0 0 ce oo 
2 
20 
2 2 3 
U U o /T "M2 * 2 | 


3. Recursive Algorithm 
The problem is solved recursively by first assuming 
the problem is solved at time K-l. Specifically it is assumed 
that the best estimate X (K-1|K-1) at time K-l and its error 
covariance matrix P(K-l|K-1) are known. 


(i) Calculate the one-step prediction 
eer (4.28 


(11) Calculate the covariance matrix for the one-step 


prediction 
T T 
P(K|K-1) = $P(K-1|K-1) 96° + GQ(K-1)G (4.29) 
(111) Calculate the prediction observation 


v(K|K-1) = HX(K/K-1) (4.30) 
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(iv) Calculate the filter gain 


G(K) = P(K|IK-1)HT([HP(K|IK-1)HT4+R(K)] | (4.31) 


(v) Calculate the new smoothed estimate 


X (K|K) X(K|K-1) + G(K) [Y(K)-Y(K|K-1)] (4.32) 


(vi) Calculate the new covariance matrix 


P(K|K) = P(K|K-1) - P(K|K-1)HT [HP(K|g-1)HT+R(K)] HP (K|K-1) 
mma  ` 
GE 
P(K|K) = [I - G(K)H] P(K|K-1) (4-33) 


In summary, starting with an estimate E 1) and its 
covariance matrix P(K-1|K-1) after receiving a new observation 
Y(K) and calculating the six quantities in the recursive 
algorithm, a new estimate X (K|K) and its covariance matrix 


P(K|K) are obtained. 


` 
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V. IMPLEMENTATION AND SIMULATION RESULTS 


In order to evaluate the three filter algorithms in a 
variety of tactical environments, an air vehicle type, and two 
tracking sensors were selected for analysis. 


P u) were selected to match 


TR D 
2 


the vehicle, and sensor statistics (on, 0%, 


Maneuver statistics (A, P 
T) were selected 
for each combination of sensor and data entry evaluated. 

One trajectory was constructed for the vehicle that consists 

of a straight track and a maneuvering track and is shown in 
Big. 39. 

The x-direction vs time of range, velocity, acceleration, 
bearing, bearing rate and bearing acceleration, were plotted 
and shown in Figs. 40-45. 

For the scenario considered, an aircraft at R = 10 KM moves 
at 100 m/sec (200 knots), can maneuver at a maximum accelera- 
tion of 4g (A = 39.24 m/sec”), and has a probability of 
maneuvering at max 0.2 (2P4 = 0.2, Pi = 0.1), and a probability 
0.5 of not maneuvering at all (P, = 0.5). Assume a lazy 
maneuver that will provide correlated acceleration inputs for 
periods between 10 and 30 sec. Hence, with an average maneuver 


I ES Er AN sec, by using Eq. (4.19), u = 1/20 = 0.05 = a, 


Il 


and the correlation coefficient P = 1 - aT L 0057. 


The two sensors were classified as: 
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A. EXAMPLE 1 - AIR SEARCH RADAR 

The radar data rate was ten samples per second (T = 10 
sec) hence, the correlation coefficient p = 1 - 0.05(10) = 0.5, 
and the sensor processing noise (measurement noise variances) 
has been taken into account, 500 m in range (or = 500 m) and 
17.4 mrad in bearing (06 = ] degree - 17.4 mrad). The 
variances of maneuvering at max acceleration, and not maneu- 
vering at all, were calculated by using Eqs. (4.17-4.18) 
4 = 


= 46193 m/sec and ao = 1.346 x 10 ` sec 


2 
ame found to be, o M2 


M1 
l. Kalman Filter Evaluation 
This model can be used with the Kalman filter 
Eqs. (4.29-4.31-4.33) to determine a set of filter gains. 
For the model assumed, by using Eqs. (4.5, 4.11-4.14, 


4.25), the following matrices were obtained: 


1 0 0 0 0 
0 T 0 0 | 0 0 
E 0 0 0.5 0 0 1 0 
0 0 È 10 0 G = 0 0 
0 0 0 1 1 0 0 
| 0 0 0 0.5 Lo 1 
1 0 0 0 0 0 
H = 
0 0 0 1 0 0 | 
Dos Tio 0 
R = SA 
0 3.02 x 10 
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34645 0 


Q = S 
0 107" 

28x10% 25x10? 0 0 0 0 

25x10” 51193 23096.5 0 0 0 

0 23096.5 46193 0 0 0 

E) = -4 -5 

0 0 0 Moo Tos 02.10 0 
E -4 -4 

0 0 0 Oo 4-10 0.673x10 

0 0 0 0 P IL ea = 


i 


E ee ed kalman Filter Evaluation 


The simplified Kalman filter is achieved by assuming 
(incorrectly) that the vehicle's change in velocity is 
uncorrelated between samples; i.e., the maneuver is white, and 
o = 0. Here, for the model assumed, by using Eqs. (4.2- 


4.4-4.5-4.26-4.27), the following matrices are obtained. 


1 10 0 0 0 0 
0 1 0 0 1 0 
Ó = G = 
0 U 1 10 0 0 
o o 5!) N 
7 4 
1 0 0 0 25x10 0 
H = R = E 

0 0 l 0 0 320210 
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A 


46193 0 


Q = Di 
0 1.346x10 
25x10* 25x10? 0 0 
25x10? 51193 0 0 
P(2|2) = _4 _5 
0 0 3.02x10 3.02x10. 
0 0 3.02x107” EE e 


~ 


B. EXAMPLE 2 - SURFACE AND AIR SEARCH RADAR 

In this example, all quantities except the following, have 
previously been defined for the alr search radar example: 
Sampling time T = l sec, hence p = l-aT = 1-0.05(1) = 0.95, 


the sensor processing noise (measurement noise variances) 


has been taken into account, Ip = 20 m, gg 7 0.1 degree = 1.74 
mrad and the variances of maneuvering were calculated and 
found to be, Gen = 461.93 m/sec”, SC = io ID Š sue 2. 


l. Kalman Filter Evaluation 


For the model assumed the following matrices were 


obtained: 
1 IL 0 0 0 0 0 0 
0 1 J 0 0 0 0 0 
0 0 0.95 0 0 0 1 0 
Ó = G == 
0 0 0 3 1 0 0 0 
0 0 0 0 L JL 0 0 
Mono o a 0 0.95 0 1 
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Le Go 400 0 


0 0 Dut 1 0 0 0 3.02x10 Š 
45.04 0 
o = lE 
0 0.13x10 
400 400 0 0 0 0 
400 1261.93 487.83 0 0 0 
0 438.83 461.93 0 0 0 
Bo) = = =6 
0 0 0 3.02x10 3.02x10 0 
0 0 0 3.02«10 ? 7.386x10 ©  1.2787x10 ? 
0 0 0 0 Te) veal UE 


2. simplified Kalman Filter Evaluation 


For the model assumed the following matrices were 


obtained: 
1 1 0 0 0 0 
0 IL 0 0 1 0 
Ó = G = 
0 0 1 1 0 0 
0 0 0 1 0 1 
1 0 0 0 400 0 
H = R = er 
oo 1º O Do 
J 
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Q = = 
0 1.346x10 
400 400 0 0 
400 1261.93 0 0 
sa) = >. Ss 
0 0 3.02x10 3.02x10 
| 0 0 2 10?  7.386=10 Š 


GE o-B FILTER EVALUATION 

The a-g filter considered in this paper, and used for 
simulation, is one of many varieties possible in this class, 
is more easily implemented than either the Kalman or Simpli- 
fied Kalman filters, and has been selected for evaluation 
since it is utilized extensively in tactical applications. 
Because it is designed to minimize the mean squared error in 
filtered position and velocity under the assumption of straight 
line target motion, it has little capability to track severely 
maneuvering vehicle. The o-8 filter examined has no provision 
to adapt to Bro ase target types, as does the Kalman filter, 
since maneuver statistics are not taken into account. The 


PAS 210-2013 -2:099—2. 38) for the a-p filter evaluated are: 


N NM N N 
Xo = Xp + o (Xy Xp) 
B; s N N 
= X == = 
Ve st pU 7 Xp) 
DW oor Red 
An = Xo + Vo T 
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where 


2(2N=1) 6 


NN: 8 N(NFT) 


anan sec or 10 sec 


in order to continually compute the least-squares line through 


the observations. 


D: COMPARISON OF FILTER ACCURACIES 

One hundred (100) Monte Carlo trials were made for each 
combination of tracking filter, tracking sensor and data entry 
procedure. 

Experimental determined filtered and predicted accuracies 
in vehicle range coordinates (range, range rate, and range 
acceleration) and bearing coordinates (bearing, bearing rate, 
and bearing acceleration) with means and variances of estima- 
tion error histories, were then calculated and plotted and 
shown in Appendix A. 

Table II shows a representative summary of the result, 
in which the prediction accuracies of each filter were 
compared on a percentage basis to that of the Kalman filter. 

The entries in the table were determined by averaging the 
experimentally obtained percentage degradations in each of 
the above coordinates. 

The Simplified Kalman filter, and the Kalman filter 
generally performed within twenty percent (20$) of each other. 
The a-ß filter performance, on the average, appears to 

be about equal to the above filters for the straight part 


of track, and about thirty to sixty percent (30-60%) worse, 
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ensor Type 






Filter Air Search Surface and 
Type Radar Air Search 
Radar 


a-B Filter 


Simplified 
Kalman Filter 


Kalman Filter 


Key: 1 = within 10-20 percent of the Kalman filter 
2 = within 30-60 percent of the Kalman filter 
Table II. Synopsis of the Accuracy Comparison 


of the Three Tracking Filters 


with the greatest degradation occurring for the maneuvering 
and accelerating part of track, because the gain vector quickly 
becomes too small to correct for the large estimation errors 


resulting from target maneuvers. 
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VI. SUMMARY AND CONCLUSIONS 


While no extensive analysis is implemented, it is con- 
Sidered that a reasonable and unequivocal comparison of the 
filters can be made from the material presented. 

The analysis of the filters and the radar system simula- 
tions presented in this paper is considered overail to be a 
Simplified but realistic model of a sophisticated system which 
could be implemented with current "state of the art" hardware. 
Based on the ensemble averages the Kalman filter obviously 
provides a somewhat better tracking response for the target 
track tested. 

The tracking ability of the o-8 and Kalman filters appears 
to be about equal for "look alike" targets in close proximity, 
under the assumption of straight line motion. The a-8 filter, 
however, provided unsatisfactory performance when the tracked 
vehicles executed maneuvers. Based also on the simulation 
results, the Simplified Kalman filter becomes attractive for 
implementation, because it provided tracking accuracies within 
ten-twenty percent of the Kalman filter. The utility of this 
filter is greatest, when either the sensitivity of tracking 
performance to assumed maneuver correlation is small, or when 
the target maneuver approaches whiteness relative to the 
sensor data rate. 

The filter implementation requirements increase in the 


following order: a-g filter, Simplified Kalman filter, 
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Kalman filter. Moreover, the "complexity factor" between the 
above filters is about two-to-one. 

Finally, in most applications, the answer to the question, 
"which filter is most accurate?", does not alone determine 
filter selection. Indeed, the following questions must all 
be answered in the filter selection process to obtain the 
"best" filter for a particular system: 

a. What are the actual accuracies of each filter? 

b. What are the relative filter accuracies? 

c. What are the tracking accuracy requirements of the 

system? 

d. How sensitive is system performance to tracking 

accuracy? 

e. What are the computer requirements of the filter? 

f. What are the computer limitations of the system? 

The list shows that filter selection involves careful balancing 
of filter accuracies, filter implementation requirements, 


and system performance goals and limitations. 
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APPENDIX A 


Simulation results, of experimental filtered and predicted 
accuracies, in vehicle range and bearing coordinates, with 
means and variances of estimation error histories, provided 
for comparison of each of the three filters. 

These results were generated using Monte Carlo simulation, 


with the following description of run: 


SE Simplified Kalman 

Filter Kalman Filter Filter 
order of system 2 d 6 
no. of measured states 2 2 2 
no. of time samples 31 31 EN 
no. of random forcing inputs 2 2 2 
no. of members in ensemble 100 100 100 
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FIG. 76. VARIANCE OF RANGE ERROR VS TIME FOR KALMAN FILTER 
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FIG. 79. VARIANCE OF BEARING ERROR VS TIME FOR KALMAN FILTER 
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